Overlapping resonances in nuclei coupling to the atomic shells 
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The resonant process of nuclear excitation by electron transition (NEET) in highly charged ions 
is investigated. In NEET, a bound electronic decay transition occurs with the simultaneous exci- 
tation of the nucleus, provided that the energies of the atomic and nuclear transition match. By 
varying the atomic charge state, the atomic transition energy can be tuned to a better match of the 
nuclear transition energy. We propose a new way to create the atomic hole in highly charged ions by 
dielectronic capture of a free electron. In order to derive the cross section for the three-step process 
composed by dielectronic capture, nuclear excitation by electron transition, and the subsequent nu- 
clear decay, a Feshbach projection operator formalism is developed. With the help of this formalism, 
the resonances in nuclei coupling to the atomic shell and the possible interference between several 
direct and resonant electronic and nuclear processes such as radiative recombination, dielectronic 
recombination and nuclear excitation by electron capture are described. 

PACS numbers: 34.80. Lx, 23.20.Nx, 23.20.-g 
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I. INTRODUCTION 

The increasing precision in atomic physics experiments 
ignited in the last decades the nuclear odyssey from the 
simplified picture of a point-like charge to its real size and 
properties. Different nuclear masses, charge distributions 
or spins have an effect on the electronic transition ener- 
gies, and the small energy corrections or splittings ob- 
served in experiments can give in turn information about 
the nucleus. Atomic spectroscopy, at present at an un- 
rivaled level of precision, can thus be used to determine 
nuclear parameters which are otherwise hardly accessible 
via nuclear physics experiments 

On the other hand, the borderline between atomic and 
nuclear physics is also the scene for several nuclear pro- 
cesses that directly involve atomic electrons. For in- 
stance, nucleus and electron can interact via the electro- 
magnetic field and undergo transitions simultaneously. A 
nucleus in an excited state that cannot decay radiatively 
transfers its energy via the radiation field to one of the 
atomic electrons which leaves the atom in the process of 
internal conversion (IC). The inverse process of IC, which 
might occur in highly-charged ions (HCI), is known as 
nuclear excitation by electron capture (NEEC) Q. In 
conjunction with HCI, atomic physics experiments in- 
volving these processes at the borderline between atomic 
and nuclear physics also open the possibility to explore 
properties of exotic nuclei 

IC and NEEC also have more exotic siblings, the bound 
internal conversion (BIC) and its inverse process NEET. 
Bound IC is a resonant nuclear decay channel which may 
occur if the nuclear excitation energy is not enough to 
ionize the atomic electron, but can induce a transition 
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between two bound electronic states. NEET is the si- 
multaneous excitation of the nucleus during an atomic 
decay transition, provided the energies of the two transi- 
tions match. There is only a limited number of nuclei in 
which such a match of the atomic and nuclear transition 
energies exist, see for instance the list of candidates for 
NEET in Ref. 0. Although it is difficult to find systems 
that fulfill the resonance condition, NEET has already 
been observed experimentally Ny], in the same year in 
which direct evidence of BIC [5| has been reported. As 
nuclear excitation mechanisms, both NEEC and NEET 
are expected to allow the study of atomic vacancy effects 
on nuclear level population and lifetimes. 

The NEET nuclear excitation mechanism has been first 
proposed by Morita back in 1973 @ as means for 235 U 
separation. This first theoretical work was followed only 
months later by a first claimed experimental observation 
[7j in 189 Os. However, these early claims of success, as 
well as the follow-up published experimental results 0- 
|Tl| were riddled with controversy. While theorists con- 
tinued arguing about the true magnitude of the NEET 
probability [12h14| , apparently overestimated in the orig- 
inal proposal Q, many studies dedicated to NEET in 
235 U in plasma environments were performed (la . Il6| . 
Other theoretical predictions confirming the important 
role that nuclear excitation mechanisms such as NEEC 
and NEET may have in dense astrophysical plasmas were 
also more recently published [13, [H| ■ 

The most recent and reliable experimental results on 
NEET have been performed on 197 Au 0, 189 0s 
[H [H| and 235 U 12. The NEET probability for the 
77.351 keV 197 Au nuclear transition that corresponds to 
an atomic K — > M transition in neutral atoms was deter- 
mined to be (5.0 ±0.6) x 10 -8 }4j. In the case of the 69.5 
keV nuclear transition in 189 Os, only an upper limit of 
4.1 x 10~ 10 for the NEET probability could be deduced 
from experiments [2l|, while for 235 IJ the occurrence of 
NEET could not be observed. 

From the theoretical point of view, a number of some- 
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times contradictory studies have been considering the 
NEET mechanism in the last decades, among which we 
mention the nonrelativistic self-consistent description for 
x-ray, Auger decay and NEET [l4|], the relativistic anal- 
ysis of NEET [13[ and a formalism developed for NEET 
in plasmas at local thermodynamical equilibrium [la ] . A 
comparison between theoretical results seeking to cl arify 
all anomalies and discrepancies is presented in Ref. }23| . 
Furthermore, the efficiency of NEET as nuclear excita- 
tion mechanism for isomer depletion has been studied 
for a number of special cases 0, 25] and was also com- 
pared or considered along with NEEC or photoexcitation 

Hi M, M M- 

The purpose of the present work is manifold. First, we 
proceed by arguing on behalf of NEET experiments with 
HCI. While atomic physics experiments with HCI contin- 
uously gained in precision, until now all the performed 
NEET experiments envisaged neutral atoms irradiated 
by synchrotron radiation light to create the necessary 
inner-shell electronic hole. Furthermore, although in the 
theoretical studies of NEET in plasmas the general im- 
portance of the ionic charge state and the environment 
has been already mentioned, we emphasize here in partic- 
ular the possibility that HCI offer to "tune" the atomic 
transition energy. A number of elements that present a 
mismatch on the order of 1 keV between the nuclear and 
atomic transition energies (E n and E a , respectively) for 
neutral atoms offer the possibility to reduce the energy 
difference (E n — E a ) by at least one order of magnitude 
for other charged states. As it will be shown in this work, 
such an improved match of the atomic and nuclear tran- 
sition energies can determine an increase of three orders 
of magnitude in the NEET probability. Considering the 
small number of energy match cases that can be found in 
nature, the partial tunability of the atomic energy brings 
additional NEET candidate elements on the list and is 
most welcome. 

Next, starting from the premises of NEET in HCI, we 
consider a new mechanism for production of the elec- 
tronic hole, namely dielectronic capture (DC). In the res- 
onant process of DC, a free electron recombines into a 
HCI with the simultaneous excitation of a bound elec- 
tron, which leaves a hole in the previously occupied 
bound shell. DC followed by x-ray decay of the excited 
electrons is known as dielectronic recombination (DR), 
an atomic physics process which is used in high-precision 
experiments with HCI to measure isotope shifts and pro- 
vide complementary data to nuclear radii obtained via 
other methods 27]. The hole created by DC can be oc- 
cupied via electronic decay that occurs simultaneously 
with the excitation of the nucleus, i.e., NEET. This two- 
step process is presented in Fig. [TJ 

By tuning the energy of the continuum electron one 
can precisely control which hole is created and which 
electronic transitions are driven. Furthermore, by using 
different types of DC, such as KLL-, KLM- or KLN- 
DC, (the notation specifies the involved atomic shells) 
one can make use of satellite electrons in higher-n shells 




FIG. 1: (1) The dielectronic capture: a free electron re- 
combines with the simultaneous excitation of a bound one. 
(2) NEET: electronic and nuclear transitions occur simulta- 
neously. The nuclear levels are depicted in the yellow panels. 

to fine-tune the atomic transition energies. As a further 
advantage, depending on the optimal charge state to min- 
imize the energy mismatch (E n — E a ), one can envisage 
a situation like in Fig. [T] in which the electronic hole can 
only be filled by an atomic transition that has the right 
energy to couple to the nuclear excitation. In contrast, 
a hole created by photoionization of an inner shell in a 
neutral atom has many decay channels corresponding to a 
multitude of electronic transitions. Furthermore, most of 
the nuclear transitions that can couple to the atomic shell 
are dipole-forbidden, and the one atomic transition which 
corresponds to NEET ought to share this multipolarity. 
In this case, the NEET atomic transition is strongly sup- 
pressed in comparison with other electric dipole {El) ra- 
diative or Auger decay channels of the atomic hole. A 
substantial gain can be thus achieved in a simple manner 
by using DC in HCI as first step previous to NEET in a 
scenario like the one presented in Fig. [T] This scenario 
may also be relevant for dense plasmas, where DR is the 
dominant electron recombination mechanism. 

We develop here a theoretical formalism that can treat 
the complex three-step process consisting of DC, NEET 
and the subsequent decay of the excited nuclear state. 
In order to derive the cross section formula for this pro- 
cess we have extended the Feshbach projector formalism 
developed and used for DR [H and NEEC @ to ac- 
count for the interaction of electronic and nuclear degrees 
of freedom. The electric and magnetic electron-nucleus 
interactions are considered explicitly. The dynamics of 
electrons is governed by the Dirac equation as required 
in the case of high-Z elements. 

Our formalism allows us to also take into account inter- 
ference processes that may occur in the sequence of res- 
onant processes. We present here a case study for 237 Np 
for which all possible interferences between radiative re- 
combination (RR) , DR, NEEC followed by gamma decay 
and the three-step NEET process under consideration are 
taken into account. The interference between the chan- 
nels of RR and DR with the sequence of DC, NEET and 
radiative decay of the excited nuclear state turns out to 
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be rather weak, but is still two to three orders of magni- 
tude larger in cross section than the NEET process itself. 
Reasons for the small interference cross sections are the 
very different time scales of the two processes, since the 
nuclear excitation needs a long time to occur as compared 
with the radiative recombination of a free electron into a 
HCI. This is not the case when considering interference 
between DC followed by NEET and NEEC, where both 
channels start from the same initial state with an elec- 
tron in the continuum and have as final state an excited 
nucleus. However, depending on the ionic configuration 
and on the recombination state for DC, this interference 
channel may not always be present. We find that also 
the magnitude of this interference term between the two 
pathways involving nuclear excitation mechanisms is the 
same as the ones involving the atomic processes RR and 
DR. At the resonance energy, NEEC turns out to be the 
dominant nuclear excitation mechanism and acts as an 
enhancement of the nuclear excitation probability. 

We would like to point out that, while for NEET 
the atomic transition matching the nuclear excitation 
is unique, NEEC has an additional degree of freedom 
for the electron capture into a bound shell as contin- 
uum electron energy and capture orbital can be varied. 
From this respect a comparison between the two pro- 
cesses needs to specify the particular conditions, for in- 
stance, the ionization degree and the capture orbital for 
NEEC. This excludes a too general argument or criterion 
for the comparison between NEEC and NEET. In the 
framework of the process sequence that we consider — 
dielectron capture in a highly charged ion followed by 
NEET and gamma decay — a unique choice for the NEEC 
channel is naturally arising when investigating the pro- 
cesses that can interfere. Previous results [l?} consider- 
ing the temperature-dependent NEEC and NEET rates 
in dense plasmas have shown that for lower ionization 
degrees, NEET is more efficient than the possible NEEC 
channels. This feature was attributed to the better over- 
lap between the bound electronic wave functions. Here, 
we show that for few-electron high-Z HCI this is not the 
case and NEEC may be the most significant contribution. 

This paper is structured as follows. We start by devel- 
oping the Feshbach formalism to account for a three step 
process consisting of DR followed by NEET and subse- 
quent nuclear decay in Section [TTJ The specific transition 
rates are identified in the cross section and their expres- 
sions are reviewed in Section Hill followed by Section ITVl 
where the interference terms are deduced and discussed. 
A numerical case study is presented in Section [V] The 
paper ends with a summary and outlook. Atomic units 
have been used unless otherwise mentioned. 



II. FESHBACH PROJECTION FORMALISM 

In this section we derive the total cross section formula 
for a three-step process consisting of DR, NEET and the 
decay of the excited nucleus. Without limiting too much 



the generality of our argument, we will treat here the ra- 
diative decay of the excited nuclear state, denoted in the 
following by 7 decay. The case of nuclear excitation fol- 
lowed by IC can be treated in a similar manner, as it was 
shown for NEEC in Ref. [29(. For simplicity we consider 
the case where the three relevant transitions energies (of 
the captured and bound electrons and of the nucleus) are 
equal. This situation is also relevant for investigating the 
interference between DR followed by NEET and NEEC. 
Furthermore, we treat the case in which the decay of the 
electronic hole occurs to the ground state, such that no 
further electronic transitions are possible. The theory 
can easily be extended to account for an arbitrary DC 
electronic configuration. 

A. Decomposition of the Fock space by means of 
projection operators 

The initial state ipi of the system consisting of the nu- 
cleus in its ground state, the free electron, the bound elec- 
trons ground state configuration and the vacuum state of 
the electromagnetic field is given by 

|^> = \I g M g ,aiJiMipm s ,0) . (1) 

Here, p is the asymptotic momentum of the free elec- 
tron and m s is its spin projection. The nuclear ground 
state is denoted by the total angular momentum I g and 
its projection M g . The bound electronic ground state 
consisting of n electrons is denoted here by the total an- 
gular momentum of the electronic configuration Ji and 
its projection Mi and any other quantum numbers on. 

In the first step DC takes place and the new state de- 
noted here as first intermediate state l^di) has the form 

\^ dl ) = \I g M g ,a 1 J 1 M 1 ,0), (2) 

where with J\ and M\ we have denoted the total elec- 
tronic angular momentum and its projection on the quan- 
tization axis, respectively, of the new electronic configu- 
ration of n + 1 electrons, and a.\ covers again any other 
quantum numbers of this state. The nucleus is still in its 
ground state and no photon has been emitted. 

In the second step NEET occurs and the system 
reaches the second intermediate state 4>d 2 of the form 

\iP d2 ) = \I e M e ,a 2 J 2 M 2 ,0). (3) 

The n + 1 electrons form now a different configuration 
with total angular momentum J 2 and corresponding pro- 
jection M 2 , and the nuclear excited state is described by 
the quantum numbers I e and M e . The excited state of 
the nucleus decays radiatively back to the ground state 
and the final state ipf is given by 

= \I g Mf, a 2 J 2 M 2 , XkLM) . (4) 

The emitted photon has the wave number fc, the to- 
tal angular momentum L and projection M. Further- 
more, A stands for electric (E) or magnetic (M) spherical 
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waves. The photonic state can be written as \XkLM) = 
a^ feLM |0), where a\ kLM is the photon creation operator. 
The corresponding conjugate annihilation operator is de- 
noted by axkLM- 

We now separate the Hilbert space into four possible 
subspaces 

• P - states with a continuum electron, 

• Qi - states with bound excited electrons and 
ground state nucleus, 

• Q2 - states with excited nuclei and ground state 
electrons (since we do not consider here the situa- 
tion when after NEET further deexcitations of the 
electronic shell may occur), 

• R - states with photon. 

We subsequently introduce operators projecting onto 
individual subspaces. Characterizing the state of the 
electron in the positive part of the continuous spectrum 
by the energy e rather than the momentum of the free 
electron, we write the projector P belonging to the first 
type of subspace as 



P = J de"£\M(M 



(5) 



For brevity we introduce the multi-index j3 to stand for 
all discrete quantum numbers of the total system. The 
projection operator of the subspace spanned by interme- 
diate states belonging to of the type |V\2i) is written as 



(6) 



with the cumulative index q\ introduced to summarize all 
discrete quantum numbers describing the excited electron 
and the ground state nucleus. The projection operator of 
the subspace spanned by intermediate states of the type 
\1pd3) is written as 



(7) 
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with the cumulative index qi introduced again to summa- 
rize all discrete quantum numbers describing the ground 
state electron and the excited nucleus. The subspace of 
the state vectors containing one transverse photon is as- 
sociated with the projection operator 



fl = E E a lkLM\ r )( r \aXkLM , 



XkLM 



(8) 



where with r we denote the quantum numbers describing 
the nucleus in the hnal state and the bound electron. As- 
suming the correction due to two or more photon states 
and due to the presence of the negative electronic con- 
tinuum are negligible, we postulate the completeness re- 
lation P + Qi + Q2 + R = 1. Similar decompositions 
involving less subspaces for the Feshbach formalism de- 
veloped for two-step NEEC and NEET can be found in 
Refs. dHH. 



B. Perturbation expansion of the transition 
operator 

The transition operator is defined as 

T(z) = V + VG(z)V , (9) 

where V is the interaction Hamiltonian discussed below 
and G(z) is the Green operator of the system given by 

-1 



G(z) = (z-HY 



(10) 



Here, z is a complex energy variable. The total Hamilto- 
nian of the system can be written as 



H — H e -\- H n -\- H r -\- H e 



(ii) 



where with H r and H n we denote the Hamiltonians de- 
scribing the radiation field and the nucleus as discussed 
in Refs. @- The subscripts n, e and r stand for the nu- 
cleus, the electrons and the radiation field, respectively. 
The electronic part is given by 28] 



H, 



N 

E 

»=i 



(ca 4 • pi + (Pi - l)mc 2 + Kuci(^)) + 



1 E-- 

(12) 

with a the vector of the Dirac matrices, pi the momen- 
tum of the ith electron and Vnuci the nuclear potential. 
Furthermore, N is the number of bound electrons in the 
ion and = \?i — fj\. We adopt the Coulomb gauge 
for the electron-nucleus interaction because it allows the 
separation of the dominant Coulomb attraction between 
the electronic and the nuclear degrees of freedom 0, 



H, 



- I d 3 r 

en 1 u 1 r 



Pn{r n ) 



(13) 



Here, p n (f n ) is the nuclear charge density and the in- 
tegration is performed over the whole nuclear volume. 
As it has been shown in Ref. [2j, this Hamiltonian term 
describes only the electron-nucleus interaction in tran- 
sitions of electric multipolarity. The interaction of the 
electron with the transverse photon field quantized in 
the volume of a sphere of radius R is given by 



H er = —3* ■ A = — 



XkLM 



(4 



XkLM 



a - A xkLM (r) + H. c 



(14) 

with the vector potential of the quantized electromag- 
netic field H 7 ! 

A(f f )= (AxkLM(r)a{ kLM + A* XkLM (r)a X kLM 

XkLM 

Here, the two independent solutions of the wave equation 
for the A X kLM(r) are 



A( M )kLM{r) 
A(E)kLM{r) 



4:7rck 



R 



3L{kr)Y? L {0, V ) 



(16) 



47Tcfc 

R 



Vx(ji(*r)lu(^)) 
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where the quantum number k is discretized by requiring 
the proper boundary conditions at a perfectly conduct- 
ing sphere of radius R. The Yff L (Q, p) denote the vector 
spherical harmonics [3lj . 

Similarly, the interaction of the nucleus with the elec- 
tromagnetic field is given by the Hamiltonian 



H nr = -- V (17) 

c 

XkLM 

x ( a XkLM J d 3 r n j n (r n ) ■ A X kLM{r n ) +H- c^J , 

where j n (r n ) is the nuclear current. 

With the help of the projector operators introduced in 
the previous subsection, we can separate the perturbation 
part in the Hamiltonian, H = Hq + V with 



H = PHP + QxHQx + Q 2 HQ 2 + RHR , 



(18) 



and 



V = H-H 

= PVQ 1 + Q 1 VP + PVQ 2 + Q 2 VP 
+ PVR + RVP + Q 1 VQ 2 + Q 2 VQ 1 
+ Q X VR + RVQ 1 + Q 2 VR + RVQ 2 . 



(19) 



The interaction Hamiltonian V is the one responsible 
for the possible transitions between the four subspaces 
described in Section III Al While H n and H r only con- 
tribute to the unperturbed Hamiltonian Hq, and H en , 
H nr and H er only contribute to the interaction V, the 
electronic Hamiltonian H e enters both expressions. The 
terms PH e P and QiH e Qi with i = 1, 2 describe the con- 
tinuum and the bound electrons, respectively. In addi- 
tion, the electron-electron interaction term in H e is also 
responsible for DC and Auger decay in the interaction 
terms PH e Qi and Q±H e P. 

The cross section of the process can be expressed by 
the transition operator as follows: 



d<jj— > r , . 2tt , , . , 



ie)\^)\ 2 Pf , (20) 



with tpf and ipi as final and initial eigenstates of Hq, re- 
spectively. This cross section is differential with respect 
to the angle of the photon emitted in the process. Fur- 
thermore, Fi denotes the flux of the incoming electrons, 
and pj the density of the final photonic states. 

We use the Lippmann-Schwinger equation to write the 
perturbation series for T(z) in powers of V with the 
Green function Gq{z) of the unperturbed Hamiltonian 
H : 



T{z) = V + VG (z)V + VG (z)VG (z)V + 



(21) 



Since the initial state of the DC process is by definition 
an eigenstate of P, and the final state after the 7 decay 



is an eigenstate of R, we consider the projection RTP of 
the transition operator: 

RTP = RVP + RVGqVP + RVGqVGqVP 

+ RVG Q VGqVGqVP + ... (22) 

Here and in the following we omit the argument z. The 
first term in Eq. (l22t accounts for RR and does not con- 
tribute to the three-step process under consideration. 
The second term RVGqVP = RH er Q 1 G n Q 1 H ee P + 
RH nr Q 2 GoQ 2 H en P accounts for the first order of two 
possible processes: DC followed by x-ray emission, which 
is DR and proceeds without any participation of the nu- 
cleus, and NEEC followed by radiative decay of the nu- 
cleus for the case of electric multipole transitions . In- 
serting the spectral resolution © of Q 2 in the second 
order in V for the NEEC term we arrive to 



y>ARrp\ii>i) = 



(ipf\H nr \q 2 )(q2\H en \ipi) 
z-E° 



(23) 



The energy E® denotes the unperturbed eigenvalue of the 
state I q). 

The third order term RVGqVGqVP can be decom- 
posed as: 



RV1Gq1V1G 1VP = 



RVPG q PVQxG QiVP 

RVPG PVQ 2 G Q 2 VP 

RVPGoPVRGqRVP 

RVQ 1 G Q 1 VQ 2 GoQ2VP 

RVQiG a QiVRG a RVP 

RVQ2G0Q2VQ&0Q1VP 

RVQ 2 G Q 2 VRG RVP . 



The sixth term 



RVQ 2 GoQ 2 VQ x G n QxVP = 

RH nr Q 2 Go Q 2 H en Qi GoQ\H ee P (24) 

7 NEET (electric) DC 

describes DC followed by NEET of electric transitions 
and contributes to the cross section of the considered 
process. Among the decomposition terms in the sum 
above we can also find the one responsible for NEEC 
for magnetic nuclear transitions 0. For DC followed by 
NEET, the term responsible for transitions of magnetic 
multipolarity will only emerge from the decomposition of 
the fourth order term in the perturbation expansion and 
has the form 

RH nr Q 2 Go Q 2 H nr RGoRH er Qi Go Q\H ee P . (25) 

7 NEET (magnetic) DC 

We have used here the result 

QH nr RG RH er P = QH^P , (26) 
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with the Hamilton operator 

H, = --a [ d\ n ^H, = -a ■ A(f) (27) 
c J \r-r n \ 

obtained in Ref. Q which accounts for a bound electron 
that decays by exchanging a virtual transverse photon 
with the nucleus. 

The fourth order decomposition in the perturbation 
expansion also contains the exchange of a virtual photon 
between a bound and a continuum photon, i.e., the Breit 
interaction in the term 




NEET {electric) DC (Breit) 



Together with the Coulomb electron-electron interaction 
term in the interaction Hamiltonian, the Breit term con- 
tributes to DC and Auger decay. Using results well 
known in the theory of DR (2|| HH we can write the 
Hamiltonian responsible for DC and Auger decay as the 
sum of the Coulomb and the Breit contributions 

1 1 \- ^ cos(wry) 
H ^ = ol^ — -l^ a i- a i — — 

a^.,-, ^ x cos(uirn) — 1 

where u) is the frequency of the exchanged virtual photon. 

Following the procedure described in Ref. |2J and 
adopting the so-called isolated resonance approximation 



we can continue with the next terms in the expansion 
up to infinite order. In Ref. 0], the infinite perturbation 
expansion was shown to introduce energy and width cor- 
rections into the energy denominator of the lowest order 
amplitude. The sum over all orders could be performed 
with the help of a geometric progression. In this case, 
since we have two projector operators Qi and Q2 corre- 
sponding to two types of intermediate states, the infinite 
perturbation expansion can be reduced to summing two 
geometrical progressions to obtain 

^ f \RT{z)P\i>i) = (30) 

Ed! d 2 (V'/ I H nr I V>rf 2 )(lpd 2 \H N \lp dl )(l/j dl \ H ee I ijjj ) 

(z~E d2 + ^T d2 )(z-E dl + ^T dl ) 

with the notation Hn = H en + H^. The energies E dl 
and E d2 of the electronic and nuclear transitions, respec- 
tively, are including the electron-electron interaction and 
the corresponding Breit corrections and all the electronic 
and nuclear self-energy contributions, as well as the nu- 
clear polarization contributions for the studied nuclear 
transition, E di = E° d . + A£f + AE% + AE*f + AE^ SE 
0. Correspondingly, the widths r rf2 (nuclear) and t dl 
(atomic) include electronic radiative and Auger widths 
for state d\ as well as nuclear radiative and IC widths for 
state g?2- 

Finally, the total cross section for our three-step pro- 
cess, after summing over all final states and averaging 
over all initial states will has form 



a(E) = 



Fi 2J, 



1 ^ ^ 

M g ,M 2 ,Mf L,M,X 



( tpf I H nr I ipdz ) ( ^d 2 1 Hn \ip dl }{4> dl \H ee \ipi 



d% ,d 2 



r 



E-E, 



di 



[E — E d 



Pi- 



rn 



C. Resonance strength for the three-step process 



The resonance strength (also known as integrated cross 
section) is obtained by integrating the total cross section 
over the whole energy spectrum. A simplification of its 
expression can be obtain by studying the relation be- 
tween the two widths T dl and r,j 2 that appear in the 
total cross section denominator. For the considered case, 
the final electronic state is the ground state. Since the 
nuclear excited state width is much smaller than the elec- 
tronic excited state width, T dl is approximately given by 
the electronic width of the initial state and T d2 is the 
nuclear width, such that T dl T d2 . Neglecting the vari- 
ation of the free electron momentum and matrix elements 
on the energy interval of interest determined by the two 
Lorentzian profiles and dominated by the very narrow 



nuclear width, the integration resumes to 

dE 



{E-E dl Y+ T -^\ ((E-E d2 )* + ^ 



2tt 



r d2 (E dl - E d2 y + -± 



(32) 



The product between the flux of the incoming electrons 
Fi and the density of the initial electronic states pi does 
not depend on the normalization of the continuum wave 
functions 32] , 



FiPi 



P 



(27T) 



(33) 
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The resonance strength of the three-step process can be 
then written as a product of three terms, 

°W dE = J? 2J i + l T d2 ' (34) 
where we have introduced the Auger decay rate 
. _ 2tt x - 

" 2JTTT ^ 

x / tm^iaiJiMtlH^laiJiMipm^pi, (35) 



the NEET probability 



1 y 

l(2Ji + l) , r ./f 



(2/ fl + l)(2Ji + l) 



(36) 



Mi,M 2 ,M c ,M g 



\{I e M e ,a 2 J 2 M 2 \H N \I g M g ,a 1 J 1 M 1 )\ 2 

( Edl -E d2 y + T -±- 

and finally the radiative decay rate of the nuclear excited 
state, 



2tt 



2L 



T £ 



M c ,M f ,\kLM 



x K/gA/^AfcLA/liJ^I/eAfe)! 2 ^. (37) 

Our formalism reproduces the results presented in 
Ref. QJ, [HI for the NEET probability for the partic- 
ular case when the second intermediate state d 2 is the 
electronic ground state. Determining the total cross sec- 
tion of the studied process requires the calculation of the 
transition rates Aa u and Ay, the NEET probability V ne et 
and the initial and final state energies. In the actual cal- 
culations we neglect the additional nuclear polarization 
corrections AE^ P to the bound electron energies [33j j . 
since their value is smaller than the nuclear transition 
energy uncertainties and not known with good accuracy. 



III. THE NEET PROBABILITY FOR 
ELECTRIC AND MAGNETIC TRANSITIONS 

The calculation of the NEET probability requires 
knowledge of the atomic levels and widths and the 
evaluation of the interaction matrix element in the 
numerator of Eq. (f36|) . In the following we sketch 
the calculation of the interaction matrix element 
{I e M e , a 2 J 2 M 2 1 Hn I IgMg, a 1J1M1) following the ap- 
proach used for the NEEC rates 

In order to calculate the NEET probability we have 
considered the matrix element of the electric and mag- 
netic interactions between the electron and the nucleus. 
We write the wave function of the system as the prod- 
uct wave function of the electronic and nuclear states. 
Unlike for NEEC, now we have in both initial and final 
state only bound electrons. The electronic multiconfig- 
uration wavefunctions \a\J\M\) are obtained by Slater 



determinants constructed with relativistic orbitals \riKm) 
in the framework of the multiconfiguration Dirac-Fock 
method (MCDF). The relativistic orbitals are eigenfunc- 
tions of the angular-momentum operators j 2 and j z and 
of the parity operator for each electron. The notation 
n accounts for the principal quantum number, k is the 
relativistic angular quantum number k = ±(j + 1/2) for 
I = j ± 1/2 such that j = \k\ — 1/2, and m is the pro- 
jection of the angular momentum operator j on the z 
axis. 

The electron-nucleus interaction Hamiltonian describ- 
ing electric multipolarity transitions can be written using 
the multipole expansion as 



E 

L=QM=-L 



Yl M (Q e ) [ 

T T J 



-j^Y L M{Qn)Pn{r n ) ■ 

The radius r e (r n ) denotes the electronic (nuclear) radial 
coordinate and f2 e (f2 n ) stands for the corresponding solid 
angle, while the Y5 m( (f2 e ) denote the spherical harmon- 
ics. We make the simplifying assumption that the elec- 
tron does not enter the nucleus, such that r e > r n . The 
Hamiltonian can then be written as 



H, 



E E 

L=0 M=-i 



4tt 



1 



2L + 1 



YZ M (Ve) — QLM (39) 



with the help of electric multipole moments Qlm = 
J d 3 r n r^YLM(^n)pn{r n )- The nuclear part of the ma- 
trix element can then be structured in the reduced elec- 
tric transition probability 



Bt(EL, I g — > I e ) = 



(40) 



Using the single active-electron approximation and fol- 
lowing the calculation in Ref. Q , we obtain for the NEET 
probability for electric transitions 



d(S) 
neet 



4tt (2J 2 + 1) 
(2L + l) 2 (2 Ji + 1) 

x c(j 2 LJi;±Q± 



Bt(EL,I g ^I e ) (41) 



{Ed x - E d , 2 



Here, C (j 2 L Ji; | |) stands for the vector coupling 



coefficients and the radial integral R L j j 2 is given by 



R 



(E) 

L,J U J 2 



drr 



-L+l 



fn2K2 ( r )/rilKl ( r ) 



9n 2 K 2 (r)9n lKl {r)] 



(42) 



where g nK (r) and f nK (r) are the large and small compo- 
nents of the bound Dirac radial wave functions describing 
the orbitals of the active electron 



l^nKm (0 
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and the indices 1 and 2 denote the two states ip^ and 
tpfei respectively. 

For magnetic multipole transitions, the magnetic 
Hamiltonian in Eq. (|27p can be written using the multi- 
pole expansion as 



H„ = -a ■ A = - 



47T 



c 2L + 1 

LM 



(43) 



Using also in this case the approximation that the elec- 
tron docs not enter the nucleus, we obtain 



LM 



4tt 



2L+ 1 



M LM d-Yll*{n 



44) 



where the magnetic multipole operator is given by 

M LM = I d 3 rr L Yll(n n ) ■ j n (f n ) . (45) 



Similar to the case of the electric transitions, all the nu- 
clear information can be contained in the reduced mag- 
netic transition probability 



Bt{ML,I g ^I e ) = 



1 



2/„ 



-\{I e \\M L \\I g )\' , (46) 



whose value can be taken from experimental data or 
from theoretical calculations considering different nuclear 
models. After some angular momentum algebra we ob- 
tain for V neet 



(M) _ 4tt ( 2J 2 + 1 W U r j _> n 



x {k x + n 2 fC ( J 2 L Ji;i ~ 



1? 



(M) 



(^ -S d2 ) 2 + ^ 

The electronic radial integral is in this case 



(47) 



R 



L,J lt J 2 



Wl = / drr- L+1 [g n2K2 (r e )f n2K1 (r) 
Jo 

+ fn 2 K 2 {r)gn 1K1 (r)] ■ 



(48) 



Given the different parity of the electric and magnetic 
multipole moments, a transition of a given multipolarity 
L is either electric or magnetic. We consider in the fol- 
lowing only the cases of transitions with a certain value 
of L which do not present mixing ratios between electric 
and magnetic multipoles of different multipolarities. 



IV. INTERFERENCE TERMS 

As already mentioned in the introduction, our three 
step process consisting of DR, NEET and 7-ray emission 



has a number of competing processes in each step that 
open a number of interference channels. Interference can 
occur whenever the initial and final states of two compet- 
ing processes are the same. For exemplification, we will 
discuss in the following a suitable case in which inter- 
ference between all electron recombination and nuclear 
excitation mechanisms may occur. This is the case for 
DC occurring in He-like HCI — in our numerical case we 
have considered capture of a free electron into He-like 
237 Np which has a suitable nuclear transition for NEET. 
The four interference pathways are presented in Figs. [5] 
andO (1) RR, (2) DC + x-ray decay, i.e., DR, (3) NEEC 
+ 7 and (4) DC + NEET + 7. Here for brevity we have 
replaced "followed" by the sign "+" . 



(1) 



2p- 
2s- 



1s- 



RR 



(2) 




DC 







2p — e— 




is c — e— 

x-ray 


• 



FIG. 2: Schematic picture of the two interference pathways 
involving purely atomic processes: the direct process of RR 
(1) and DC + x-ray decay (2). 



(3) 



2s- 

1s- 



NEEC 







2p 


2s *- 

y-ray 


— — 
• 



(4) 







'-t— 1 

?p — 1 — e- 




DC 





2s- 
1s- 



NEET 







2P 


2s 1 

7-ray 


— e- 
• 



FIG. 3: Schematic picture of the two interference pathways 
that couple the atomic and nuclear degrees of freedom: the 
two step process NEEC + 7 decay (3) and the three-step 
process of DC + NEET + 7 decay (4). 

In the common initial state ipi we have a He-like ion 
in the electronic ground state Is 2 with its nucleus in the 
ground state, and a free electron of energy Eg and an- 
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gular momentum quantum numbers of the partial wave 
\enjm). The final state consists of the ground state Li- 
like ion, ground state nucleus, and an El photon of en- 
ergy matching the nuclear and atomic transitions. As 
long as the state ls2s2p 3 / 2 decays directly to the ground 
state (consistent with the ls2s2p 3 / 2 —> ls 2 2p 3 / 2 transi- 
tion rate being several orders of magnitude smaller than 
the ls2s2p 3 / 2 — > ls 2 2s one), the four pathways cannot 
be discerned and thus quantum interference occurs. 

The Hamiltonian terms that describe the four interfer- 
ence pathways can be found in the perturbation expan- 
sion of the transition operator presented in Section III Bl 
Except for the direct process of RR, which is described by 
the term RVP, all the other three combinations include 
resonant processes that can only be accounted for in the 
perturbation expansion by summing the corresponding 
terms up to infinite order. This procedure, originally de- 
veloped for the calculation of DR cross sections, has been 
also adapted for NEEC and for NEET in this work. 

Thus, we can write the transition operator of the four in- 
terference pathways that connect the initial and the final 
states ipi and tpf as 

(il3 t \BT{z)P\il3i) = (4> f \RH er P\^) 

(ip f | H nr | ip d > 2 }(<p d , 2 \H N \ tp^ ) | H ee | ifjj ) 



£ 

di dL 



(z-E d , 2 + ±r d , 2 )(z-E d[ +±T d{ ) 



+ £ 

d 2 

+£ 



(tpf | H nr 1 1pd 2 ) ( lpd 2 | H N | tpi 



z-Ea 



2 L d 2 



(ij} f \H er \^ dl )(ij: dl \H ee \ipi 



z- E, 



i T dl 



(49) 



By introducing the expression above in the cross section 
formula in Eq. (|20p. we obtain a sum of ten terms, four 
of them corresponding to the cross sections of the four 
considered processes pathways and the other six to the 
interference between them. In the following we label the 
simple terms with and the interference terms between 
pathways i and j as <xy, with i,j = 1,2, 3, 4. We are now 



interested in the three interference terms that involve the 
process under consideration, labeled as channel (4) . Note 
that two of the other three interference terms, <j\ 2 and 
<7i3, have been already investigated in Refs. [28ll35|. 

In the following we will address the test case of 237 Np 
and sketch the calculation of the three interference terms. 
The NEET nuclear transition is occurring between the 
ground state and the excited state at 102.959 keV. As it 
will be shown in the numerical results section, this energy 
is best matched by the El atomic transitions between the 
ls2s2p 3 / 2 state and the ground state ls 2 2s of the Li-like 
ion. The active electron is undergoing a transition from 
the 2p 3 / 2 to the Is orbital. The decay of the 2s electron 
to the if -shell occurs much slower and thus justifies the 
considered one active electron scenario. The DC occurs 
in the ground state He-like ion, such that the free electron 
is captured in the 2s orbital with the simultaneous exci- 
tation of one of the Is 2 electrons to the 2p 3 / 2 state. Since 
this is an El transition, both RR and DR processes have 
large cross sections, such that the corresponding interfer- 
ence term in the cross section is expected to be relevant 
for the total cross section. The numerical study in the 
next section shows to what extend this is the case. 



A. RR interference term 

In order for interference between RR and pathway (4) 
to occur, the continuum electron that recombines via RR 
or DC has to have the same initial and final energies 
and quantum numbers. In both cases, the recombination 
occurs into the 2s orbital. In RR the final electronic 
bound ground state is reached (the nucleus remains in 
the ground state all the time), while for the three-step 
process (4) the recombination into the 2s orbital leads to 
the excited state ls2s2p 3 / 2 that decays via NEET. 

By introducing the expression of the transition opera- 
tor matrix element in Eq. (1491) in the expression of the 
total cross section, we obtain the interference cross sec- 
tion term between pathways (1) and (4), 



a u (E) = ^PiPf E T \ £ 27TT £ 2~TTl £ f {a2 J2M2 ' ( E ) kLM \ H er\mJiM iP m s )* (50) 

P M,Mf,M 2 m s 1 Mi 9 M g \ 

(I g Mf, {E)kLM\H nr \I e M l e ){I e M' e ,a 2 J 2 M 2 \H en \I g M g ,a' 1 J[M[)(^^^ 

X k*&M> (E-E d ,+ir d ,)(E-E d2 + ir d2 ) 



Here we have assumed that NEET can occur only via a 
specific nuclear level denoted by the total angular mo- 
mentum I e (whereas its projection is not fixed and thus 
must be summed over). Furthermore, for the sake of 
brevity, the non-participating electronic or nuclear wave- 



functions in the matrix elements of 7 decay or DC, re- 
spectively, were omitted in the notation. The partial 
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wave expansion of the continuum wave function reads 
\pm s ) = i^YLS^) 



x C [I - j;mi m s 



\EKjm) 



(51) 



where e is the energy of the continuum electron measured 
from the ionization threshold, e = \Jp 2 c 2 + c 4 — c 2 . The 
orbital angular momentum of the partial wave is denoted 
by I and the corresponding magnetic quantum number 
by mi, while the partial wave phases A K are chosen so 
that the continuum wave function fulfills the boundary 
conditions of an incoming plane wave and an outgoing 
spherical wave. The total angular momentum quantum 
number of the partial wave is j — \k\ — 5. The ex pre ssion 
of the RR matrix element (see, for instance, Refs. [35ll36j ) 
is given by 



(a 2 J2M2, {E)kLM\H er \eKjm) 
J 1 



Anck 
R 



C(j L J 2 ;m M m d ) 



2 j + 1 

4"7T 



2 



R 



RR 

J2,j,L ■ 



(52) 



The Wigner 3j-symbol used in the expression above is 
related to the Clebsch-Gordan coefficients via 



C(ji 32 j;mi m 2 m) = (-l)™+h~h ^2j + l 



3i 32 3 
mi m,2 —m 



(53) 



The expression Rj^j l involving the radial integrals is in 
this case given by 



R 



RR 

J2,j,L 




Y (L + 1)(2L+1) 
with the radial integrals 



(LI^-^-k)!^) (54) 



+ (*d - k)I+ +1 ) 



It 



drr 2 j L (kr) (g„ 2K2 (r)/ eK (r) ± g £K (r)/„ 2K2 (r)) , 

(55) 

containing the spherical Bessel functions jL(kr) and the 
large and small radial components of the relativistic con- 
tinuum electron partial wave function g eK (r) and /^(r), 
respectively. 

For the matrix element of the interaction Hamiltonian 
between the nucleus and the radiation field H nr , we make 
the assumption that the wavelength of the radiation is 
large compared to the nuclear radius, kR <C 1, so that 
the Bessel functions that appear in the expression of the 
field vector potential (|16p can be approximated in the 
first order in kr as 

(kr) L 



j L (kr) 



(56) 



(2L + 1)!! 

In this case the electric solution of the wave equation can 
be written as 



A 



(E)kLM 



(f) = - 



47rcfc ^/(L + 1)(2L + 1) 



(kr) 



R 

-lyM 
1 LL 



(2L + 1)!! 



(57) 



With the use of the continuity equation for the nuclear 
current j n we obtain for the matrix element 

(I g M f ,(E)kLM\H nr \I e M e ) 



(-1) 



Z.-Me + l 



VL + T 



Aitck 

R 
ik L 



C (I g I e L; M g - M e M) 



V^(2L+1) (2L+1)!! X 9 



IgWQLWh 



(58) 



Finally, one can relate in a simple manner the DC matrix 
element to the Auger rate for the case under considera- 
tion, making the observation that since the total angular 
momentum of the bound Auger final state Is 2 is 0, there 
is only one corresponding partial wave for the recombin- 
ing continuum electron with 3 = 1/2 (j = 3/2) and k = 1 
(k = —2) for the J\ = 1/2 (Ji = 3/2) intermediate state. 
We then have 

(axJxM^He^aiJiMijm}^^) = 
(aiJiWHeeWlaiJifij^S^j^M^i , (59) 

and the Auger rate is given by 

A Au = 2ir\(a 1 Ji\\H ee \\[a l J i j}j 1 )\ 2 p l . (60) 

Combining the expressions above, we obtain for the in- 
terference cross section the expression 



<Tli(E) 



16tT 



Pi B\(E\,I g ->/ e )(2J 2 + l)(-l) 



k L+1 / L + l ( J 2 j L 
(2L + l)!!yL((2L + l) 3 ^ | -| 

R< l E Lj 2 R Tj,l («i J il I H ee 1 1 [a, j a ] . h ) (61) 
1 \ 



x Re 



(e — E dl + §r rfl ) (e — E d2 + §r d2 ) 



B. DR interference term 

DC is the first step for both pathways (2) and (4). 
By DC, the doubly-excited electronic state ls2s2p 3 / 2 is 
reached, that can subsequently decay either by NEET 
or by x-ray emission. NEET may occur only when the 
2p electron decays to the Is hole, while the 2s electron 
acts as spectator. In principle, the 2s electron may also 
decay by x-ray emission to the isf-shell and contribute 
to DR, but this process is much slower. We therefore 
consider in the following the interference between NEET 
followed by 7 decay and El x-ray emission of the doubly- 
excited electronic state following DC. We can write the 
interference term in the cross section as 
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, , (2tt) 4 1 T„l^ 1 \ 1 

1 M,M f ,M 2 "ia Mi 9 M g 

( Eg, , h Ml (<*2 J 2 M 2 , (g) kLM\H er I oi J! Mi ) * (ai Ji Mi |ff ee I a, J l M l ym s )* 

X I j 

Y — Ed 1 — 2^di 

(I g Mf, {E)kLM\H nr \I e M' e ) (I & M' e , a 2 J 2 M 2 \H en \I g M g , a[ J(M()(ai J[M[ \H ee \ai JiMifim s 



(62) 



*E E 



[e - + §r d /j (i? — E d . 2 + §r d2 ) 



+ H. c. 



Using the expressions of the four matrix elements given in Eqs. (j4Tj) . (|52|) . (|58p . (|59|) . and the summation properties 
of the Clebsch-Gordan coefficients, we arrive to the interference cross section term 



16tt 3 
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<J2i(E) = ^ Pl Bt{ElJ g -> 7 e )(2J 2 + l)(2j + l)(-l) 3/ <= 



h/„ fc L+1 I L+l fj 2 L j 
(2L + 1)\\V L(2L+lf y -| 



4f i , J 2 ^ J 2 ",Zl I (« 1 J l 1 1 ^ee 1 1 [<*i Jj] J, ) | 2 RC 



1 

[E — Ed 1 — |r dl ) (£7 — £? dl + §r dl ) (E — Ed 2 + §r d3 J 

I 



.(63) 



The radial matrix element i?j 2 L corresponding to 
x-ray emission has the same expression as Rj 2 R j i in 



Eq. (|54l) with the difference that both initial and final 
electronic states are bound. The continuum electron par- 
tial wave functions g ere (r) and f eK {r) are therefore re- 
placed by g niKl (r) and /„ lKl (r), respectively. 



cleus. Instead of DC, the free electron can recombine 
into the bound shell with the simultaneous excitation of 
the nucleus in the process of NEEC. The excited nuclear 
state decays in both channels radiatively. Thus, the in- 
terference occurs between DC followed by NEET on the 
one hand and NEEC on the other hand. 



C. NEEC interference term 

Finally, the last interference term of interest is the one 
that involves in both channels the excitation of the nu- 



J 



The interference term in the cross section is given by 



aa j E ) = ^p tPf Y — [ dn - y y y 

34V 1 p 2 F F) ^ 4tt J p 2 ^ 2Ji + 1 ^ 21 Q + 1 ^ 



E 



M,M f ,M 2 m 3 ' Mi M g 

(I g M f ,(E)kLM\H nr \I e M e )* (I e M e ,a 2 J 2 M 2 \H en \ I g M g , a, J, M,pm s ) * 



e — Ed 2 — f r d2 



E E 

M' B a\J{M[ 



(I g M f , [E)kLM\H nr \I e M' e ) (I e M' e , a 2 J 2 M 2 \H en \I g M g , a l 1 J[M[)^ 1 J' 1 M[\H ee \a l J l M l V™s 



(e — Ed 1 + §r dl ) [E — Ed 2 + §r d2 



+ H. 



(64) 



NEEC and NEET are described by the same Hamiltonian 
H en . However, the matrix element for NEEC denotes in 
the electronic part a transition from a continuum state 
to a bound state, i.e., recombination, while for NEET 
we have a transition between bound states. Using the 
continuum electron wave function expansion in partial 



waves and after performing the summations we obtain 



cr 3i (E) 



64tt 4 



Pi Bt(El,I g -> I e )Bl(El,I e -> I g ) 



X [(2L + 1)!!]2 6 \ J2 L J > 2 ° 2 ) R ^,J, R l M 

L + l (aiJiWHeeWlaiJij]^) 
X L(2L + ir {E -E d2 y + r -¥ 

x(2J 2 + l)Re' ' 



(E — E^ + iTd, ) 



(65) 
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NUMERICAL RESULTS 



related through the formula 



In the following we present in detail a numerical ex- 
ample for the process of NEET in HCI. One of the inter- 
esting candidates is 237 Np with its 102.959 keV El nu- 
clear transition which comes close to the L3 — > K atomic 
transition. The mismatch between the atomic and nu- 
clear transition energies is listed in Ref. |3j as 1.89 keV 
for neutral atoms. This energy difference enters the de- 
nominator of the NEET probability in the second power, 
as showed in Eq. (|36|) . Consequently, a reduction of the 
mismatch can enhance substantially the NEET proba- 
bility, up to the maximum of AE = in which case the 
denominator in Eq. (|36[) will be determined by the width 
of the first intermediate state F^ , i.e. the width of the 
electronic state. 

The first step is therefore to identify the optimal elec- 
tronic configuration for NEET and to calculate the tran- 
sition energies between the 2p and Is subshells for differ- 
ent charge states and electronic configurations. For the 
calculation of bound electron level e nerg ies and widths 
we have used the GRASP92 package 37]. The acronym 
GRASP stands for General-purpose Relativistic Atomic 
Structure Program and is a suite of FORTRAN codes 
for various calculations of relativistic atomic structure. 
The MCDF approximation is used for the calculation 
of atomic stationary states and transitions among them 
38]. GRASP92 is an improvement of the previous ver- 
sions and includes approximate QED corrections for the 
electronic energy levels. The finite size of the nucleus, 
i.e., its radius Ro, is also considered in the GRASP92 
wave functions and has a sensitive effect on the energy 
levels of the bound inner-shell electron. The accuracy of 
the GRASP92 bound electron energies is expected to be 
on the level of less than 10 eV for few-electron systems 
such as He-like and Li-like ions, well within the exper- 
imental uncertainties. The calculated atomic transition 
energies and widths are listed in Table U for a number of 
ion configurations. We see that the energy mismatch is 
smallest for the Li-like configurations. 

Choosing the most advantageous configurations as the 
[Is2s2p 3 / 2 }j set, we proceed to calculate the NEET prob- 
ability and the resonance strength for the three-step pro- 
cess of DC followed by NEET and 7-ray emission. The re- 
sults, including continuum electron energies, Auger rates, 
the NEET probability and the resonance strength are 
presented in Table [TTJ The Auger rates were obtained 
with a computer code extension of GRASP developed 
for DR calculations [39l - l4l| . The radial wave functions 
that enter the integral in Eq. (|42]) are obtained with the 
GRASP92 code. The nuclear radiative rate was calcu- 
lated according to the formula [3(| 

A ^ E1) = lZl + iw E2 T Bi{XL ' h -+ /g) ' (66) 

with L = 1. The two reduced transition probabilities for 
the emission, respectively the absorption of a 7 ray are 



Bl(\L,I e ^I g ) = ^ft±B J t(\L,I g ^I e ) . (67) 



91 2I e + 1 
The width of the excited nuclear state is then 



E 



(68) 



where Aic is the IC rate for the decay of state d 2 and v is 
the sum over all decay channels. In our case, due to the 
high binding potential of the two X-shell electrons, only 
IC of the 2s electron is possible. The internal conversion 
coefficient for this state is small, 2.14 x 10~ 2 42| for the 
102.959 keV transition, such that the nuclear width is de- 
termined by the radiative decay of the nucleus. It should 
be noted that the 102.959 keV nuclear state does not only 
decay directly to the ground state, but also via two inter- 
mediate states at E = 59.5409 keV and E = 33.196 keV, 
such that both radiative and IC decay rates sum up three 
contributions v. The most important contribution for the 
nuclear width is given by the radiative and IC decays to 
the 59.5409 keV state. 

The calculated NEET probabilities are small, on the 
order of 10 -9 . Compared to previous relativistic calcula- 
tions that considered NEET occurring in neutral atoms 
irradiated by x-rays to generate an inner-shell hole [l3j], 
we see an enhancement of more than three orders of mag- 
nitude of the NEET probability. This is due to the very 
small energy mismatch AE = E n — E a for the case of 
NEET in HCI. After taking into account the different 
energy mismatch AE — E n — E a in the denominator of 
'P ne et considered in the previous works, we find that our 
values for the NEET matrix element confirm the results 
in Ref. [ill [23| and are three orders of magnitude smaller 
than the results in Ref. 14]. 

The cross section as a function of the free electron en- 
ergy has the shape of a narrow Lorentzian profile, with 
the width determined by the excited nuclear state width 
Td 2 = 2.75 x 10~ 6 eV and the peak value at the nuclear 
transition energy E =102.959 keV of a = 1.8 x 10~ 4 b 
for the [ls 1 2s 1 2pi /2 ] 3/2 and a = 6.2 x 10~ 5 b for 

[ls 1 2s 1 2py 2 ]i/2 initial configurations, respectively. 

After investigating the magnitude of the NEET cross 
section, we turn now to the competing processes denoted 
in Section [W] as channels (1), (2) and (3) and the in- 
terference terms between them. For 237 Np, the nuclear 
transition is electric dipole allowed, such that the com- 
peting channels of RR (direct recombination into the 2s 
orbital) and DR (DC in the [ls 1 2s 1 2pg , 2 ] j state followed 
by x-ray decay) are strong. The total cross sections for 
channels (1) and (3) were calculated following the ap- 
proaches presented in Refs. @, H, |36|. For the NEEC, 
RR and x-ray rates we have evaluated the radial inte- 
grals i?g j2 , R* R jL and i?™ L that enter Eqs. (HI, 
(|6"3")l and (|6"5"|) numerically. Relativistic Coulomb-Dirac 
wave functions for the continuum electron and GRASP92 
bound electron wave functions were used. 
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TABLE I: The electronic transition energy E a , mismatch to the nuclear transition energy AE — E n — E a and excited electronic 
state width T c i 1 for various ion configurations as initial state for NEET in 237 Np. 

Configuration E a (keV) AE (keV) F dl (eV) 

[ls 1 2s 2 2pj /2 2p^ /2 ] 1/2+ 101.664 1.295 70 

[ls 1 2s 2 2p? /2 2p| /2 ]i_ 101.701 1.258 76 

[ls 1 2s 1 2pi /2 ] 1/2 _ 102.879 0.080 30 

[ls 1 2s 1 2pJ /2 ] 3 / 2 _ 102.999 -0.040 22 

[ls 1 2pl /2 ] 1 - 103.279 -0.320 29 



TABLE II: NEET probabilities and total resonance strengths for the three-step process of DC into the [Is2s2p 3 / 2 ]j doubly- 
excited state followed by NEET and 7 decay of the nuclear excited state. The atomic transition energy E a , continuum electron 
energy Ep and intermediate state widths are given, together with the Auger rate Aa u - 



J 


E a (keV) 


Eg (keV) 


r dl (eV) 


r d2 (eV) 


Aau (1/s) 




Ja{E)dE (beV) 


1/2 


102.879 


69.353 


30 


2.75xl0~ 6 


3.36xl0 12 


7.2xl0~ 10 


8.19xl0 -8 


3/2 


102.999 


69.275 


22 


2.75xl0~ 6 


2.52xl0 12 


2.8xl0~ 9 


2.36xl0~ 7 



Due to a rather small nuclear reduced transition proba- 
bility B(E1) 53, the NEEC resonance strength for elec- 
tronic capture into the 2s orbital with the simultaneous 
excitation of the 102.959 keV nuclear state is not very 
high, only S = 6.21 x 10~ 4 beV. Just as for the case of 
RR, in the calculation of the NEEC rate we have consid- 
ered only the recombination of the continuum electron 
with the partial wave (k = 1 for the configuration with 
J = 1/2 and k — —2 for the configuration with J = 3/2, 
respectively) that contributes to the interference. The 
NEEC resonance strength is still three to four orders of 
magnitude larger than the NEET resonance strengths 
presented in Table HU This can be traced back to the 
overlap between the bound and continuum radial wave 
functions for Np 91+ . We have investigated and compared 
the magnitude of the radial integrals R L ■ j2 and R L ' Jx j2 
for several values of the continuum electron energies. It 
turns out that for HCI, where the continuum electron 
recombines into the strong field of the ion, the radial in- 
tegrals involving the continuum electron wavefunctions 
are larger than the ones involving only bound electronic 
states. Among NEET bound-bound radial integrals, we 
find that the overlap is largest for 2s — > Is transitions. 

The energy dependence of the cross sections for the 
four recombination channels is decisive for the behav- 
ior of the interference terms. Both channels (3) and (4), 
which involve the excitation of the nucleus, present a very 
narrow Lorentzian energy profile centered on the nuclear 
transition energy. Although the resonance strength, i.e., 
the integrated total cross section is very small, the nar- 
row width of the Lorentzian profile (determined by the 
nuclear excited state width) makes the cross section at 
resonance to be comparatively large: 1CP 4 and 10~ 5 b 
for the three-step process involving NEET and 143 b for 
the two-step process involving NEEC. 

The atomic recombination processes, on the other 



hand, do not present such narrow resonances. RR is a di- 
rect process with no resonance profile and for the energy 
interval of interest around the nuclear transition energy 
the cross section can be approximated as constant at the 
value 27.6 b. For the total RR cross section, two par- 
tial waves in the continuum electron expansion contribute 
to the cross section. However, only one of these partial 
waves can participate in the interference. The contribu- 
tions of the two partial waves in the total cross section 
are 6.83 b for k = 1 (capture into [ls 1 2s 1 2pgy 2 ] 1 / 2 ) and 
21.8 b for k = —2 (capture into [ls 1 2s 1 2py 2 ] 3 / 2 )- 
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FIG. 4: Interference terms in the cross section an (RR), 
024 (DR) and (734 (NEEC) as a function of the continuum 
electron energy for recombination into Np 91+ . See text for 
further explanations. 

DR is a resonant process whose width of the Lorentzian 
profile is determined by the width of the excited elec- 
tronic state, on the order of tens of eV. Consequently the 



14 




E (keV) 

FIG. 5: Total cross section and (74(E) (NEET) for the recom- 
bination of a free electron into Np 91+ at the resonance energy 
for DC and NEEC. See text for further explanations. 



resonance cross section is only 2.44 b and 2.38 b for the 
[ls 1 2s 1 2py 2 ] 3 / 2 and [ls 1 2s 1 2p\^ 2 ]i/ 2 configurations, re- 
spectively, at the continuum electron energy correspond- 
ing to the nuclear resonance, and stays constant on the 
whole narrow energy region of interest. Thus, although 
the total resonance strength of both RR and DR is or- 
ders of magnitude larger than the ones of the processes 
involving nuclear excitation, their values at the nuclear 
resonance energies are not very large. For the case of 
237 Np, the RR and DR cross sections at the resonance 
energy are orders of magnitude larger than the one of 
NEET, but much smaller than the one for NEEC. The 
interference terms au(E), (124(E) and 034 all inherit 
the very narrow width of channel (4) , and the magnitude 
of the interference cross section is determined by the reso- 
nance cross sections of the individual processes. The RR 
and DR contributions in the interference are expected to 
be largest for electric dipole allowed transitions, as it is 
the case of 237 Np. In Fig. |4]we plot the three interference 
terms considering for DC the [ls 1 2s 1 2py 2 ] 3 / 2 _ configu- 
ration as a function of the continuum electron energy 
on a very narrow interval around the resonance energy. 
We see that both the purely atomic processes as well as 
NEEC have interference terms with similar peak values 
at the resonance energy, on the order of 10~ 2 b. This 
value, although small, is two to three orders of magni- 
tude larger than the NEET cross section value. We note 
that this situation is qualitatively different than the one 
presented in Ref. (3|| for the interference between RR and 
NEEC. In that case, compared to the large peak value of 
the NEEC cross section, the interference term was orders 
of magnitude smaller and therefore negligible. 

The total cross section of process (4) with and with- 
out including the interference terms and the competing 
processes are presented in Fig. [3] for DC leading to the 
[ls 1 2s 1 2py 2 } 3 / 2 _ state. The cross section plateau is de- 
termined by the RR and DR cross sections, while the 



peak value at the resonance energy is determined by 
NEEC. Surprisingly, although RR and DR have much 
larger resonance strengths than the processes involving 
nuclear transitions, it is NEEC which helps rise the to- 
tal cross section of the three-step NEET process at the 
resonance energy. Due to the very large resonance cross 
section of channel (3), we see that among the interference 
terms, the term 034 (E) is smaller but comparable with 
<7i4 (i?) and (724(E). However, more than the interference 
cross section term, it is the NEEC term itself that is 
the relevant contribution for nuclear excitation in the re- 
combination process. We conclude that in HCI, the two 
nuclear excitation mechanisms NEET and NEEC, when 
both possible, might deliver contributions of similar im- 
portance and should be treated on equal footing. 



VI. SUMMARY 

In this work we have investigated a new aspect of elec- 
tron recombination into HCI involving the coupling of 
the atomic shell to the nucleus in the process of NEET. 
Our scenario involves the resonant process of DC into 
HCI to create the electronic hole needed for NEET. HCI 
present the advantage that the atomic level energies are 
very sensitive to the ion charge state and offer the possi- 
bility to optimize the match between atomic and nuclear 
transition energies. The NEET probability can thus be 
enhanced by several orders of magnitude compared to 
neutral atoms. This can be of great advantage for exper- 
iments aiming at investigating NEET. Furthermore, both 
HCI and DC are predominant in dense astrophysical plas- 
mas where NEET is expected to play an important role 
for the population of nuclear excited states. 

We have developed a versatile formalism for describ- 
ing complex processes actively involving atomic electrons 
and nuclei and derived total cross sections and resonance 
strengths for the direct and resonant channels of elec- 
tron recombination. The total and interference cross sec- 
tion terms for the processes of RR, DR, DC followed by 
NEET and 7 decay and NEEC followed by 7 decay were 
deduced and their magnitude investigated for a test case. 
Our results show that for HCI NEEC may be the most 
important contribution to the total recombination cross 
sections and that the interference terms, although small, 
are still larger than the NEET cross section. This is qual- 
itatively different from case of the typical NEET scenario 
investigated until now in experiments, where the specific 
inner shell vacancy is created in neutral gold atoms 197 Au 
4] and NEEC cannot occur simultaneously with NEET. 
The impact of our findings for nuclei of astrophysical rel- 
evance in isomer depletion and nucleosynthesis is still to 
be investigated. Calculations for such particular cases as 
well as estimates about the feasibility of a new type of 
NEET experiment with HCI are in progress. 
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